PROPERTIES OF NEARBY INTERSTELLAR HYDROGEN DEDUCED FROM 

LYMAN a SKY BACKGROUND MEASUREMENTS G. E. Thomas 

ABSTRACT For a sufficiently rapid relative motion of the solar system and the nearby interstellar gas, 
neutral atoms may be expected to penetrate the heliosphere before becoming ionized. 
Recent satellite measurements of the Lyman a emission above the geocorona indicate 
such an “interstellar wind” of neutral hydrogen emerging from the direction of 
Sagittarius and reaching to within a few astronomical units of the sun. We present here a 
detailed model of the scattering of solar Lyman a from the spatial distribution of neutral 
hydrogen in interplanetary space. This asymmetric distribution is established by solar 
wind and solar ultraviolet ionization processes along the trajectories of the incoming 
hydrogen atoms. The values of the interstellar density, the relative velocity, and the gas 
temperature are adjusted to agree with the Lyman a measurements. The results may be 
interpreted in terms of two models, the “cold” model and the “hot” model of the 
interstellar gas, depending on whether galactic Lyman a emission is present at its maxi- 
mum allowable value or negligibly small. The cold model predicts densities between 0.03 
and 0.06 cm" 3 , temperatures below a few hundred degrees K, and a wind velocity of 
about 10 km/sec. The hot model predicts densities between 0.09 and 0.12 cm" 3 , tempera- 
tures in the range 10 3 to 10 4 °K, and a wind velocity of about 6 km/sec. Since it is 
unlikely that any appreciable galactic emission is present, the results suggest that the solar 
system is within a large turbulent cell of hot tenuous gas, the so-called intercloud 
medium. 

INTRODUCTION 

We are concerned with the small part of the energy 
emitted from the quiet sun that is dissipated in the local 
interstellar medium. This energy is mainly in the form of 
extreme ultraviolet radiation (EUV) below the Lyman 
continuum (k < 911 A) and the fast -moving protons in 
the solar wind. The deposition of both types of energy 
results in ionization of the local medium. The interstellar 
neutral hydrogen, whose abundance is of the order of 
90 percent, undergoes photoionization and charge 
exchange with solar wind protons. The spatial distribu- 
tion of the resulting ionization about the sun is a sensi- 
tive function of the relative velocity of the sun and the 
local gas. For the special case of zero relative motion, an 
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extensive H II region will be formed similar to the classi- 
cal Stromgren sphere. The radius of this region depends 
on the interstellar hydrogen density rtjj and has been 
calculated by Williams [1965] to be ~1500«|_j“‘ 2/3 in 
AU. A balance between ionizations and recombinations 
is established on a time scale of the order of 10 6 years 
[Brandt, 1964]. 

For such a large time scale, it is a simple matter to 
show that a relative velocity of as little as 500 cm/sec 
can destroy this equilibrium state. For a relative velocity 
of the order of 1 km/sec or greater, the distribution of 
ionization may be described by a dynamic equilibrium 
[Blum and Fahr ; 1910a ] . The basis of this description is 
the probability of ionization P(r) of an interstellar 
hydrogen atom at the distance r from the sun and is given 

dsJ(r')IVJr 0J (1) 
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where J(r ) and V r (r) are the total ionization rate and 
the relative velocity at the point r . The integration is 
performed over the particle trajectory to the point r. 
The neutral density is calculated by adding up the con- 
tributions from all trajectories passing through the point. 
Detailed calculations have been performed by Blum and 
Fahr [1970a], hereafter referred to as BF, who show 
that for a relative velocity of 20 km/sec an asymmetrical 
ionization cavity is formed. The small dimension of the 
cavity lies in the forward direction of flow and is of the 
order of 1 to 5 AU, depending on solar activity. 

The BF theory does not include the effects of the 
bow-shock termination of the solar wind, which prob- 
ably occurs well outside the effective boundary of the 
ionization cavity [Ax ford etal , 1963; Hundhausen, 
1968; Blum and Fahr ; 1969] . Blum and Fahr show that 
neutral hydrogen flows through this region with negligi- 
ble interaction. The countereffect of the ionization 
cavity on the location of the shock front has been 
recently studied by Semar [1970] . His results suggest 
that this interaction is not likely to be important for the 
range of interstellar densities determined in this paper 
(0.03-0.12 cm' 3 ). Accordingly, these regions may be 
considered as separate and independent and we will 
refer hereafter only to the ionization cavity. 

How can such an ionization cavity be observed? It 
was first suggested by Kurt [1967] that the Lyman a 
sky background emission represents the scattering of 
solar Lyman a from beyond the cavity. Kurt and 
Germogenova \ [\961] calculated the diffuse radiation 
field within a spherical cavity surrounding the sun as a 
function of the cavity radius. Using the measured inten- 
sity of Lyman a emission from the Venus probe Zond 1 
of 75 R(\R = 10 6 /47t photons cm"“ 2 -sec~ 1 ) they obtained 
a value of 3500 AU. 

The more recent Lyman a measurements of 
Chambers etal [1970] have been cited as evidence by 
Blum and Fahr [1970&] for their model of a cavity of 
much smaller dimensions. These measurements showed a 
smoothly varying Lyman a emission with a maximum in 
the direction nearest the apex of the solar motion. 
(Hereafter we refer to this direction as simply the apex.) 
It is located at the coordinates RA = 270°, 5 = 30° 
[Alien, 1964] . A minimum was observed at a direction 
about 180° from the maximum. As pointed out by Blum 
and Fahr this is the variation of scattered solar Lyman a 
expected from hydrogen outside a cavity elongated in 
the direction of the apex. The inflowing hydrogen atoms 
would approach most closely in the direction of motion. 
They would be most likely to be ionized in the opposite 
(downwind) direction. 


We have recently reported extensive Lyman a back- 
ground measurements taken from a spinning OGO satel- 
lite [Thomas and Krassa , 1971; Bertaux and Blamont, 
1971] . Our intensity amps clearly show that the 
Lyman a intensity reaches its greatest value in the 
direction RA as 263°, 8 =~22°. (It is now clear that the 
measurements of Chambers etal [1970] could not 
define its precise location but only a meridian of maxi- 
mum brightness, since all their scans were made over a 
single great circle in the sky.) This suggests that the 
nearby hydrogen possesses its own peculiar velocity, 
which combined with the solar system apex motion 
produces the observed direction of the “interstellar 
wind.” 

When we compared two sets of measurements 
separated by an interval of 7 months, we discovered a 
parallax effect in the intensity distribution (recently 
confirmed from a preliminary examination of recent 
OGO data taken in September 1970). This finding is the 
most direct evidence of a nearby source of scattering and 
makes it possible to triangulate to determine the 
effective distance of the scattering region. 

We will describe a theory of the distribution of 
hydrogen in the solar system, similar to that of Blum 
and Fahr, but including two important effects neglected 
by them. The first of these is the outward solar Lyman a 
radiation pressure force, which nearly balances the 
inward gravitational force [Wilson, 1960 \ Brandt, 1961], 
causing the inflowing hydrogen atoms to pass the solar 
vicinity in straight-line trajectories. Furthermore, the 
velocities of the hydrogen atoms at all solar distances are 
unvarying and equal to the relative velocity of approach 
at infinity. This results in a considerable mathematical 
simplification of the theory. It also disposes of the 
necessity to account for the absence in all the measure- 
ments of a cusp-like feature predicted by BF as a result 
of gravitational focusing in the downwind direction. 

The second effect we consider is that of a random 
(thermal) component of velocity in the interstellar 
medium. In the absence of gravitational focusing, the 
OGO measurements demand a considerable “filling-in” 
of the downwind direction. If this is not provided by an 
external source of diffuse emission from the galaxy, then 
a rather high temperature (up to 10 4 °K) of the nearby 
medium is needed. As discussed in the final section, such 
a high temperature is in agreement with 21-cm measure- 
ments of the temperature of the low-density intercloud 
medium. 

THE OGO DATA AMD THE INFERRED MOTION 
OF THE LOCAL HYDROGEN 

Figures 1 and 2 are maps of Lyman a isointensity con- 
tours for the periods September 12-13, 1969 (SU 1) and 
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Figure 1 LAST contour map of the Lot sky background in celestial coordinates for SU 1. The units labelling the 
contour lines are Rayleighs. The light curves indicate the paths of the field of view. The path of the earth between 
2000 GMT \ 12 September 1969 , and 1000 GMT, 13 September 1969 , as seen from the satellite is shown as a dotted 
line. The data used for constructing the map were all taken in this period. The dashed lines are the contour lines 
which have been interpolated beneath the geocorona. 
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Figure 2 Contour map of the Lot sky background in celestial coordinates for SU 3. The path of 
the earth is shown between 0600 GMT, 1 April 1970, to 2200 GMT, 2 April 1970 . 




April 1-2, 1970 (SU 3). These data were taken by the 
Laboratory for Atmospheric and Space Physics at the 
University of Colorado (LASP). The Paris maps [Bertaux 
and Blamont , 1971] are very similar. Figure 3 is a com- 
parison of the two sets of data plotted along the ecliptic 



Figure 3 The LASP and Paris Lyman a data plotted 
along the ecliptic I equator for the two spin-up periods 
SU 1 and SU 3. 


equator. This plane passes very nearly through the points 
of maximum and minimum Lyman a intensity and is a 
convenient coordinate plane for displaying the impor- 
tant spatial variations. We do not consider the overall 
disagreement of about 30 percent significant in view of 
absolute Calibration difficulties in the ultraviolet. The 
discrepancies in the location and shape of the maximum 
and minimum are probably a result of minor problems in 
the removal of unresolved stars from the background in 
the LASP data, and in the removal of effects of charged 
particles and statistical noise in the Paris data. 

The important features of the data for this discussion 
are the shifts along the ecliptic in the positions of the 
maximum and the minimum that occurred between the 
two dates. For the LASP data these displacements 
amount to 37° and 20°, implying effective distances of 
3 AU and 6 AU for the source regions of the maximum 
and minimum, respectively. The Paris data show a larger 
displacement of the maximum (50°) but indicate a 
motion of only 10° for the minimum that is opposite to 
that expected for a parallax motion. Because of the 
greater difficulty of accurately locating the position of 
the minimum, it is possible that this effect is not real. 
Indeed, there is a 13° discrepancy between the two sets 
of data for the SU 3 period in locating the maximum. 
This suggests that errors of at least this magnitude must 
be present in the low brightness regions of the maps. 


Accordingly, we will place more emphasis on the 
parallactic motion of the maximum. 

The mean position of the Lyman a maximum is 
located at RA = 263°. This is assumed to be the true 
direction of the “interstellar wind” V r , which is the 
result of two velocities Vjj and-V 0 shown in figure 4; 
Vjj is the velocity of the interstellar gas and V 0 is the 
velocity of the solar system (20 km/sec), both relative to 
the local standard of rest (LSR). The wind speed Ffj and 
the angular separation | of the vectors V r and Vfj are 
given by 


v H =V e 2 + V r 2 - 2V rV Q cos? (2) 


V r -V Q cos? 

cos?= ^ (3) 


| is the angle (53°) between the mean position of the 
Lyman a maxima and the apex. In a later section we 
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Figure 4 The geometry of the interstellar wind. The 
relative velocity V r is the difference between the hydro- 
gen velocity Vjj and the solar system velocity V Q . The 
ecliptic plane passes through and is perpendicular 
(within 10° ) to the line connecting the Lyman oc maxi- 
mum and} minimum ( the direction of V r ,). 

show that V r is d'f the order of 5 to 10 km/ sec. From 
equation (2) the value of Fjj is calculated to be about 
19 km/sec for both values of V n with corresponding 
values for f of 111° and 96°. This indicates that the 
“true” velocity of the interstellar hydrogen is directed 
from the vicinity of the constellation Triangulum 
Australe, making an angle of about 30° with the galactic 
plane. 


671 




At first sight this conclusion appears to contradict the 
findings of Venugopal and Shuter [1967] and previous 
radio investigations. On the basis of 21-cm line shifts 
from nearby hydrogen they found no systematic differ- 
ential motion between neutral hydrogen and the LSR. 
However, the average distance of the observed hydrogen 
emission can be estimated to be MOO parsecs from the 
galactic latitude of their measurements (±20 to ±45°). 
Even though small on a galactic scale this distance 
(~10 7 AU) is still enormous compared with the size of 
the region (~10 AU) probed by solar Lyman a scatter- 
ing. It is possible that 21-cm measurements represent an 
average over a large number of irregularities in the 
velocity field. 



THE DISTRIBUTION OF INTERPLANETARY 
HYDROGEN 

We will refer to as “interplanetary” those nearby inter- 
stellar hydrogen atoms that have survived ionization 
and make an appreciable contribution to the scattering 
of solar Lyman a in the vicinity of the earth. The 
Lyman a scattering plays two roles: it renders the 
neutral hydrogen visible to ultraviolet photometers and 
it produces a force that effectively balances the gravita- 
tional force. To demonstrate the latter role, we note that 
the average force on a hydrogen atom is the product of 
the photon momentum hvjc and the number of photons 
scattered per second. The latter quantity is known as the g 
factor and is useful in expressions relating an airglow 
intensity to the column density of scattering atoms 
along the line of sight [< Chamberlain , 1961]. Barth 
[1969] calculates this value to be 2.1XHT 3 sec* 1 . We 
have used the high-resolution solar Lyman a measure- 
ments of Bruner and Reuse [1969] and the photon flux 
reported by Hinteregger [1970] to calculate a slightly 
larger value g = 2.3X1 0" 3 sec -1 . For the latter value, the 
outward radiation force and the gravitational force differ 
by only 2 percent. 

The "Cold" Model of Interstellar Hydrogen 

We first consider the case of a unidirectional streaming 
velocity, the “cold” model. With no gravitational bend- 
ing of the trajectories, the continuity equation for the 
density of neutral hydrogen is 


Figure 5 The coordinate system for describing the 
interplanetary distribution of hydrogen . The X Y plane 
is the ecliptic plane . The earth’s positions at the two 
spin-up periods SU 1 and SU 3 are shown. The inset 
shows the addition of the thermal velocity Vj to the 
relative velocity of flow (i.e the mean flow) 

J= f k -h F X dX (5) 

k ^ is the ionization cross section of atomic hydrogen 
and is the solar photon flux. In equation (4) a term 
of order V r /V sw i-s ignored. Here ^sw is the (constant) 
solar wind velocity (^300 km/sec). As we will show 
V r < 30 km/sec and hence the ignored term is 
<10 percent. 

Assuming that and F \ vary as r" 2 over the region of 

interest, equation (4) may be integrated along the 
streaming direction from +°° to the point r (r,d 7 X). 
The result is 



d (a4> + J)n H 

£ n ^ ~V r (4) 

where x is the coordinate measured along the streaming 
axis (fig. 5), o is the charge exchange cross section, and 
4> is the solar proton flux. The quantity J is the 
photoionization rate per atom and is given by 


% and J e are, respectively, the solar wind flux and 
photoionization rate at r e = 1 AU; w j|(°°) is the inter- 
stellar hydrogen density; r, 0, and X are the sun-centered 
coordinates of the point r shown in figure 5 . Since the 
measured positions of the maximum and minimum are 
nearly in the ecliptic plane we will identify y with the 
ecliptic longitude measured from the meridian passing 
through the +X axis (ecliptic longitude, 263°). Similarly, 
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the angle 5 shown in figure 5 is the ecliptic latitude 
measured from the earth. For points in the ecliptic 
(XY) plane (5 = 0), equation ( 6 ) reduces to expressions 
given by Semar [1970] and Tinsley [1971] . 


Scattering of Solar Lyman a by 
Interplanetary Hydrogen 

Given the distribution of hydrogen the 

single-scattered Lyman a intensity is calculated from 


where 


n H 


= n H (°°)exp 




r c f(d)~ 

r 


( 9 ) 


m 


e_ 

sin 0 


(0 < 0 < tt/2) 


7T 0 — 7t/2 
2 sin 0 


(tt/2 < 0 < 7f) ( 10 ) 


The quantity r c is the “effective cavity radius” in the 
forward direction (+JY) of streaming. It depends on the 
wind velocity and the total loss rate but is independent 
of rt H (°°). 


4nl = ge dsn u r e 2 /r 2 (14) 

J o 

where ds is an element of length along the direction in 
which the intensity is to be calculated, and g e is the 
g factor at the radial distance r e (1 AU)> For the direc- 
tion outward from the sun (7 = a, 8 = 0 ) the integration 
over s may be performed analytically for both models: 

“Cold” model 


4t rl(r 0 ) = 


y e 2 ”H<°°) 

r c f(Q) 


1 -exp 



(15) 


The "Hot" Model of Interstellar 
Hydrogen 

We now include in the theory the possibility of a ran- 
dom component of velocity superimposed on the 
streaming velocity. It is sufficient to assume that in the 
moving frame a single velocity Vj is distributed isotropi- 
cally. Since the mean free path of hydrogen is ~10 3 AU 
the atoms will describe collisionless trajectories. For 
simplicity the equations will apply to points in the XY 
plane. The generalization to three dimensions is straight- 
forward. The density at any point is due to contribu- 
tions from atoms whose velocities lie within a cone of 
half-angle 0 , where 


(V T IV r ) sin 0 
tan * ” ~ l-(pyF,.)cos/3 


( 11 ) 


where j 3 is the angle between the velocities V r and V^ 
shown in figure 5. The resulting V r + Vy has a com- 
ponent along the streaming direction of 


V x = (V r 2 -2V r V T cos 0+ V T 2 ) 1 12 (12) 


Since the random velocity Vj’xs, distributed isotropically 
each stream (defined by 0 ) is weighted only by its proba- 
bility of ionization. The total is therefore a sum over all 
angles j3, that is 


« H (°°) 


Z H ' 


27T 


/; 


d \ 3 exp j 


l-r c v r f[e-mv - 

rV x 0) 


( 13 ) 


“Hot” model 


4 nl(r 0 ) = 


Ve 2w H(°°) 

2nr c 



1 - exp[-^((3)//- 0 ] 

m 


where 


m= 


v r f[e-m\ 

(16) 

v x 

(17) 


For more general directions of viewing, it is necessary to 
perform a numerical integration of equation (14). 


COMPARISON WITH OGO LYMAN a. 
MEASUREMENTS 

The adjustable parameters of the theory are ftp^ 00 ), 
Vj iV n and r c . The important observational quantities 
are the values of the Lyman a intensity, the ratio of the 
maximum and minimum intensities, and the parallactic 
displacements of the maximum derived from measure- 
ments taken at different periods. The uncertainties 
attached to each of these quantities (insofar as the LASP 
data are concerned) are estimated to be: (\) intensity— a 
factor of 2 is not unlikely in the laboratory calibration 
[Pearce et ai , 1971]; (2) maximum/ minimum ratio — 
estimated to be accurate to within 5 percent, this ratio 
involves the relative accuracies, together with the small 
degree to which the effect of stars influences the deter- 
mination of the sky background; and (3 ) parallactic dis- 
placement - a maximum or a minimum is defined by a 
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smooth curve through at least three points on three dif- 
ferent circular scans of the sky (as in fig. 3). The maxi- 
mum error in this quantity is of the order ±5°, causing 
the error in the parallax measurement to be about ±7°. 
If similar errors are present in the Paris data, it is clear 
the discrepancies of figure 3 are not unexpected. Rather 
than carry the estimated formal uncertainties through 
each of the derived quantities of the theory, we prefer to 
be conservative and assume that these quantities will lie 
somewhere in the range defined by the LASP and Paris 
data. 

The procedure in fitting the cold model to the OGO 
measurements is to vary the parameter r c in the intensity 
calculations that apply to the periods SU 1 and SU 3. 
The ecliptic longitude of the earth on those dates was 
351° and 191°, respectively, corresponding to values of 
a (see fig. 5) of— 88° and +72°. From the intensities 
calculated on the two dates, the angular separation of 
the intensity maxima is obtained and is plotted in 



I 5 <0 15 

r c {AU)-~ 

Figure 6 The separation in degrees along the ecliptic of 
the Lyman a maxima for the periods SU 1 and SU 3, 
denoted by “parallax” versus the parameter r Q1 the 
“ effective cavity radius. ” The horizontal lines show the 
measured values of the separation for the two OGO 
experiments. 


figure 6 as a function of r c . The value implied for the 
LASP data is r c = 6.4±1.0 AU, and for the Paris data, 
r c = 5.1 AU, with a similar uncertainty. 

Figure 7 is a plot of the density for the cold model in 
the 0=0° direction as a function of r. This also shows 
the “source function”77j^>)/r 2 , which is the integrand in 
equation (14) for the direction a = 0°, y = 0°. This illus- 
trates that the effective distance of the emitting direc- 
tion is about 3 AU, the value obtained from a simple 
geometrical parallax calculation. 

For the direction a. = 1 80° , 7=180°, the intensity 
calculated from the cold model should be zero, a conse- 
quence of complete ionization for this direction. This 
can be seen in equation (9), since /(0)->°°as’0 -^ 7 r.For 
the actual positions of the earth (and the OGO satellite) 
the predicted intensity is finite but extremely small for 
this general “downwind” direction. To account for the 
measurement of 240 R in this direction, it is necessary 
to invoke an outside Lyman a source. The most likely 
candidate is the diffuse Lyman a radiation from the 
galaxy. Adams [1971] has considered the detailed his- 
tory of Lyman a photons originating in recombination 
in H II regions with subsequent scattering in H/ regions 
and the possibility of absorption on interstellar dust 
grains. His upper limit of a few hundred Rayleigh for the 
diffuse sky background corresponds to the case of no 
dust absorption. Therefore, although it is unlikely, we 
cannot absolutely exclude 240 R of isotropic galactic 
Lyman a as a possible contribution to the OGO measure- 
ments. When this value is added to the prediction of the 
cold model and the total normalized to the LASP 
measurement at the maximum, the value for /W 00 ) is 
0.061 crif 3 . 

For the hot model, the intensity in the downwind 
direction is provided by scattering from neutrals that 
enter the “shadow” region by virtue of their inclined 
trajectories. For Vj> V r , a wind from all directions 
occurs. For this case, neutral atoms may approach the 
solar system from the opposite direction of the flow. 
For V t » V r , the wind flux and thus the scattered 
intensity is nearly isotropic. The observed 
maximum/minimum ratio of 2.38 (for the LASP data) is 
achieved with a thermal velocity equal to 1 .44 times the 
wind velocity. With the quantity VjjVy known, the par- 
allax may now be calculated from equations (13) and 
(14) as a function of r c . The results are shown in 
figure 6, from which the values of r c may be deduced as 

before. Finally, the parameter flj^ 00 ) is adjusted until the 
maximum intensity is matched. The interplanetary den- 
sity for the hot model is plotted in figure 8 for both the 
upwind (0 = 0°) and downwind (0 = 180°) directions, 
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Figure 7 The distribution of neutral hydrogen n^(r)and the source function for Lyman a scatter- 
ing, n^rj/r 2 , plotted against radial distance from the sun for the direction of maximum Lyman a 
scattering . 



Figure 8 The distribution of neutral hydrogen for the 
“cold” and “hot” models as a function of radial distance 
from the sun . The symbols + and 0 indicate positions 
where the interstellar density is reduced by the amounts 
50 percent and ( 1/e ), respectively. 


and is compared with that of the cold model. Figure 9 is 
a comparison of the LASP measurements and the theory 
for both models. 

With the knowledge of the ionization rates and J , 

and the quantity r c , we can determine the wind speed V r 
from equation (7). Furthermore if V r is known, the ran- 
dom velocity Vj may be determined from the ratio 
Vj jV r . Using the data of Hinteregger [1970] for the 
ionizing flux F ^ between 91 lA and 280A, and the cross 
section k \ given by Allen [1964] , the value of J e is 
calculated to be 6.89X 10~ 8 sec" 1 , corresponding to a 
10.7 cm solar radio flux of 144 (in units of 
1 0‘ 2 2 watts nf 2 Hz'* 1 ). The average value for F l0 .i for 
the 6 months preceding each of the two spinup periods 
was 148 [Solar-Geophysical Data , 1970] ; thus, no scal- 

ing to a different radio flux is needed. (A 6-month 
average was selected since the interplanetary atoms will 
move through about 2 AU during that time, the approxi- 
mate width of the source function region shown in 
figure 7.) 

Typical quiet-time solar wind fluxes of 
1 .6X 1 0 8 cm" 2 sec" 1 are reported by Hundhausen [1970] . 
Increases in flux as large as 6X 1 0 8 cm" 2 sec -1 can occur 
during the passage of a disturbed magnetic sector region 
[Neugebauer and Snyder , 1966] . However, these tran- 
sient effects are expected to exert very little influence 
on the average ionization rate. The charge-exchange 
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ECLIPTIC LONGITUDE (DEGREES) 


Figure 9 The comparison with the LASP data of the theoretical Lyman a scattering intensity for 
the “cold” and “hot” models of the interstellar gas. The cold model contains an isotropic intensity 
of 240 R. 


cross section is about 2X 1 (T 1 5 cm 2 [FiteetaL, 1962]. 
Table 1 is a summary of the final results using 
a4> = 3.2X 1CT 7 sec" 1 . The temperatures for both sets of 
data are derived from the relative flow velocities V r and 
the ratio V j fV r - 1.44, using the formula 
m H P T 2 is the mass of the hydrogen atom 

and k is Boltzmann’s constant. The reader may scale 
these results according to any total ionization rate he 
chooses. 

ADDITIONAL EFFECTS NOT CONSIDERED IN 
PRESENT ANALYSIS 
Temporal Effects 

A perfect balance of radiation pressure and gravitation is 
not expected. Furthermore the net force will vary during 
the solar rotation period in response to 27-day changes 
in the line-center region of solar Lyman a \ Meier, 


Table 1 Basic heliosphere parameters for “hot” and 
“cold” models. 


Experiment 


Deduced 

quantity 


Cold model Hot model 


r c (AU) 

6.4 

10.9 

V r (km/sec) 

9.1 

5,3 

r(°K) 

0.0 

3.56X10 3 

« H (cm' 3 ) 

0.061 

0.12 

r c (AU) 

5.0 

8.5 

V r (km/sec) 

11.6 

6.9 

r(°K) 

0.0 

5.87X10 3 

( cnf3 ) 

0.030 

0.090 
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1969] . The net force may also vary appreciably along 
the trajectory due to a Doppler shift of the resonant 
frequency within the self-reversed solar Lyman a line 
profile. The use of a constant net radial force (assumed 
here to be zero) and constant values of the EUV and 
solar wind fluxes are justified because of the small dis- 
tances covered by the interstellar hydrogen atoms 
(0. 1-0.2 AU) during one solar rotation period. This pro- 
duces a strong averaging effec t on the photo ionization as 
discussed by Blum and Fahr [1970c] . 

Secondary Ionization 

In addition to the loss processes considered here, neutral 
hydrogen will be lost by impact ionization from fast 
protons and neutrals (the latter created by charge 
exchange). The measured value of the relevant cross sec- 
tion for 9 keV hydrogen atoms is 6X 1 QT 1 7 cm 2 [Allison, 
1958] , only 3 percent of the charge-exchange cross sec- 
tion. This process is therefore of negligible importance. 

Doppler Shifts 

For relative velocities ranging from 5 to lOkm/sec, the 
center of the absorption profile will be shifted by as 
much as 0.02 to 0.04A, respectively. (The earth’s orbital 
velocity of 30 km/sec can add up to ±0.1 2A to the 
measured line shift.) The high-resolution solar Lyman oc 
line profile of Bruner and Reuse [1969] shows that such 
shifts can increase the resonance g factor over its 
line-center value by only a small amount. This could 
cause a slight sharpening of the intensity maximum pre- 
dicted for a constant g factor, a desirable gain for the 
theory in view of the measured sharpness shown in 
figure d. However, it would also lead to a broadening of * 
the minimum, probably to a greater extent than the data 
indicate, although more measurements are needed in this 
portion of the sky before this can be stated with cer- 
tainty. Such modifications would be most important for 
the cold model and for the higher values of V r . 

Additional Solar Wind Loss Processes 

Hundhausen etal. [1968] have suggested that the asym- 
metrical charge-exchange process 

He ++ + H-> He + + H + 

can modify the ionization state of the solar wind by 
conversion of outflowing He ++ to He + . The measured 
cross section for this process is ^3X 1 Cf 1 6 cm 2 [Fite 
etal, 1962]. Even with the admixture of as much as 
15 percent He ++ in the solar wind [Hundhausen et ai. 


1967] , the loss rate of neutral hydrogen by this 
mechanism would not exceed 3 percent of the 
charge-exchange loss. 

Hundhausen and his coworkers also suggested that 
the anomalously high concentration of He + measured in 
the solar wind could be explained by this process if the 
incoming neutral hydrogen flux were ~2X 1 0 6 cm' 2 sec -1 . 
However, none of our models gives fluxes greater than 
IX 10 s cm" 2 sec -1 even at large distances from the sun, 

Photoionization of H 2 , OH, and H 2 0 in 

the Interstellar Medium 

Although the concentrations of these molecules are 
expected to be very small in the interstellar medium, it is 
of some interest to estimate the effects of the breakup 
of these molecules in the inflowing gas. Consider first 
the effects of EUV and solar wind on an H 2 molecule. 
Since charge exchange occurs readily and photoioniza- 
tion is more likely than photo dissociation, an H 2 
molecule will be converted mainly into an H 2 + ion 
rather than into two hydrogen atoms. These ions will be 
accelerated very quickly by the pXB electric field and 
swept away by the solar wind. ThusH 2 would probably 
not contribute to the interplanetary hydrogen atom 
population, unless it has a much greater abundance than 
H in the nearby interstellar medium . 

The H 2 0 molecule is readily dissociated in the solar 
ultraviolet with a photodissociation rate of about 
l(T s sec _1 [Anderson, 1971] . Thus, the breakup of water 
molecules occurs at distances ten times greater than 
hydrogen atom ionization. The interplanetary distribu- 
tion of H atoms would be nearly the same as if the 
interstellar medium were originally composed of H 
atoms and OH radicals. 

The ultraviolet absorption cross sections of the OH 
molecule have not been measured. However, if the inter- 
planetary hydrogen were primarily a dissociation 
product of OH (or equivalently of H 2 0), interplanetary O 
atoms would be present in comparable numbers to H 
atoms. This possibility could be tested by searching for 
the 13 04 A sky background emission. Unfortunately, the 
g factor is only 1.5X KT 6 [Barth, 1969] due to a very 
faint solar 1304A line. The long-wavelength channel of 
the LASP ultraviolet photometer [Thomas and Krassa, 
1971] could have detected an intensity greater than 
20 R. This imposes only the weak restriction that 
the interplanetary atomic oxygen is less than about 
1 cm" 3 . 

The above discussion is intended only to show that 
the measurements themselves do not necessarily exclude 
photodissociation of interstellar molecules as the origin 
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of the observed interplanetary hydrogen. To definitely 
rule out such possibilities it is necessary to calculate a 
detailed spatial distribution of the “hot atoms” 
produced by photodissociation and the Lyman a scat- 
tering expected from such a distribution. This would 
permit the assignment of upper limits to the abundances 
of hydrogen-bearing molecules in the nearby gas. 

Solar Lyman a Attenuation 

At a sufficiently great radial distance, the solar Lyman a 
will be appreciably depleted in parts of the line so that 
solar radiation pressure no longer plays an important 
role. Therefore, the incoming atoms will initially be 
accelerated. This acceleration continues until the atom 
approaches sufficiently close so that the attenuation is 
no longer effective, or until its Doppler shift causes the 
resonance wavelength to move into a spectral region | 
where attenuation is unimportant. The resulting distribu- 
tion of velocities with distance will cause an originally 
cold gas to scatter the solar line as if it were at a finite 
“effective temperature” T e . The value of T e may be 
estimated by first calculating the distance R x at which 
the line-center optical depth is unity for a variety of gas 
temperatures. The distance R\(T) ranges from 8 AU at 
T= 10 2 °K to 20 AU at T- 10 4 °K. We next assume 
that the incoming atoms are accelerated by solar 
gravitation beyond Rx (and maintain a constant 
speed and direction for r < R x ). The velocity 
increments A V (R x ) are then computed from 
( V Y 2 + 2M q GJRi) U 2 V r , where V r is the initial rela- 
tive velocity, M 0 is the solar mass, and G is the gravita- 
tional constant. The velocity increment AV may be 
converted to an effective temperature T e = m^AV 2 /2k 
as a function of R x . At the point where T e (R x ) = T(R X ), 
we obtain the velocity increment and effective temper- 
ature that is reached by the accelerating gas before it 
becomes balanced by solar radiation pressure. For an 
initial velocity of 10 km/sec (see table 1), AV is calcu- 
lated to be ~3 km/sec, and T e ~ 600° K. For an initially 
finite gas temperature the effect is to add AV to the 
thermal velocity. This argument is valid for the 0=0° 
direction only. For 0°, the velocity increment along 
the radial direction and hence T e will be somewhat less. 

The major conclusion that was drawn from the cold 
model was that the unidirectional streaming velocity 
could not fill in the downwind direction with neutral 
atoms. This is not quite true with the above modifica- 
tions, since those trajectories for 6^0° will be slightly 
bent toward the sun and will produce a nonzero popula- 
tion along the 0 = 180° axis. However, our calculations 
show that unless V T >V r (T>5X 10 3 °K for 
V r = 10 km/ sec) the effect of filling in the shadow is 


very small as far as the Lyman a scattering intensity is 
concerned. 

Multiple Scattering of Lyman a Photons 

Since the mean free path for scattering of Lyman a 
(8-20 AU) is not much greater than the effective 
dimension of the source region (fig. 7), the effects of 
multiple scattering may not be negligible. A rule of 
thumb for estimating its importance at a distance r is to 
compare the number of solar photons N s available for 
scattering to the number Ng available for rescattering 
from the surrounding gas. If N s » Ng, multiple scatter- 
ing may be ignored. At 1 AU, N s = l - 4X10 10 
photons cm" 2 sec" 1 and the value of the mean scattered 
intensity Ng ^ 4.5X 1 0 8 photons cm” 2 sec" 1 ; N s decreases 
as f 2 out to the attenuation distance R x and then 
rapidly falls to zero. According to equation (15) or (16) 
the outward first-order scattered intensity initially falls 
rapidly with r, and then approaches a variation of r l . 
However, we may estimate the average (over all solid 
angles) fails off approximately as r~ l at all distances. The 
value of r for which N s =lQNg is therefore about 
4 AU. When N s = 2Ng, r ^ 20 AU. Thus, in the source 
region of the Lyman a scattering (1-8 AU), second-order 
scattering is likely to be only 10 to 20 percent of the 
primary scattering. At greater distances, the effect 
becomes more important, but the number of photons 
available for scattering is very small. Thus we do not 
consider this effect a serious one in our analysis. 

DISCUSSION 

If the galactic Lyman a emission is as high as 200 to 
240#, the interstellar hydrogen density is between 0.03 
and 0.06 cm" 3 . In this model, the interstellar gas tem- 
perature is low (~10 2 °K) and the interplanetary hydro- 
gen is completely ionized in the downwind direction. On 
the other hand, if the galactic contribution is negligible, 
the gas temperature must be high (up to 10 4 °K) to fill 
in the downwind region with neutral hydrogen. An inter- 
mediate value of the galactic emission leads to a “warm” 
model of the gas with 0.03 < < 0.12 cm" 3 . Knowl- 

edge of the galactic emission is not very important in 
determining the interstellar density ; however, it is criti- 
cal in determining the gas temperature. Adams [1971] 
has suggested a means of differentiating the galactic and 
local components by measuring the line width of the 
Lyman a sky background. The galactic component will 
be broad (1 to 3 A). The present analysis predicts a width 
of the interplanetary component to be of the order of 
0.05 A. Its Doppler shift (up to 0.1 2A) depends on the 
direction of viewing and. the season of observation. 
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Previous determinations of the interstellar H/ density 
from radio measurements give average values of 0.4 to 
0.7 cnf 3 [see Brandt etal., 1971 for a recent discussion 
of these measurements] . Jenkins [1970] has calculated 
an upper limit of 0.04 cm" 3 for the average hydrogen 
density along the 100-pc distance to ol Vir from rocket 
measurements of the Lyman a interstellar absorption by 
A. Smith. A rocket-borne spectrometer measurement of 
an enhancement of 1216A emission from a Bootes 
(Arcturus) has been recently reported by Rottman et al 
[1971] . Assuming the emission line width of Lyman a is 
comparable to that of the sun (~lA), we can calculate 
the column density of hydrogen that would barely 
extinguish the emission line. Over the 1 1-pc distance, the 
average hydrogen density would be no larger than 
0.05 cm' 3 . For distances of a few hundred pc where a 
number of bright ultraviolet stars are available for 
Lyman a absorption measurements, the average densities 
are on the order of 0.1 cm' 3 [Jenkins, 1970] . The 
well-publicized discrepancy between the 21 -cm results 
and the Lyman a absorption results has apparently not 
yet been resolved. Our results lend support to the 
Lyman a absorption measurements; however, it should 
be emphasized that our value refers to a very much 
smaller distance scale than either of the above tech- 
niques. Small-scale irregularities in density could easily 
account for the discrepancy between our result and the 
radio measurements. Irregularities are apparently incap- 
able of resolving the problem of the 21 -cm versus the 
Lyman a absorption measurements [ Jenkins , 1970] . 

The 21 -cm observational evidence for the thermal 
state of the intercloud medium indicates temperatures 
between 1X10 3 and 5X 10 3 °K [Rohlfs etal , 1969] . 

A current theory of the heat balance of H7 regions 
predicts two stable phases of the neutral gas existing in 
pressure equilibrium [Spitzer, 1968; Field et al, 1969] . 
These ideas can account for the known existence of 
dense, cold clouds (i n 10 cnf 3 , T~ 1G 2 °K) and a 
tenuous, hot intercloud medium 0.1 cm" 3 , 

1 0 3 — 1 0 4 °K) if the medium is postulated to be 
heated by a large subcosmic ray flux. Another theory 
favors the mechanism of heating by soft X rays [Werner 
etal , 1970]. A time-dependent model involving the 
possible contributions from supernova outbursts has 
been advanced by Bottcher et al [1970] . 

Since the hydrogen clouds occupy only 10 percent of 
the galactic volume, we would expect from pure chance 
that the solar system is contained within the intercloud 
medium. Our results for the density of the nearby gas 
clearly indicate that this is the case. Our results are also 
consistent with the intercloud gas temperatures deduced 


from radio measurements, provided the galactic Lyman 
a contribution is not too large. 

SUMMARY 

Our recent Lyman a sky background measurements have 
provided much new information on the nature and loca- 
tion of the scattering sources. They indicate that most if 
not all the emission is interplanetary in origin, and place 
the 1 effective scattering distance within about 3 AU of 
the sun. Through the use of a simple model of the dis- 
tribution of the neutral interplanetary hydrogen we have 
shown that the OGO data are consistent with an inter- 
stellar hydrogen density between 0.03 and 0.12 cnf 3 , an 
interstellar gas temperature as high as 10 4 °K, and an 
interstellar wind velocity between 5 and 10 km/sec. The 
interstellar hydrogen velocity relative to the local stan- 
dard of rest is calculated to be about 19 km/sec, directed 
from the southern hemisphere at an angle of about 30° 
with the galactic plane. 

The deduced motion of the nearby gas may be a 
manifestation of a turbulent motion on a scale that is 
probably much greater than a mean free path (10 3 AU) 
but much smaller than the thickness of the galactic plane 
(10 7 AU). The gas temperature required of the “hot” 
model is consistent with 21 -cm evidence of the high 
temperature of the galactic intercloud medium. However 
the measurements are also consistent with a cold model 
of the gas, provided about 200 R of galactic emission is 
present in the sky background. Removal of this uncer- 
tainty will require high spectral resolution measurements 
above the geocorona, or measurements made from space- 
craft traveling to the outer boundaries of the solar 
system. 
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T. E. Holier I have two comments. The first is with regard to the Lyman a that you DISCUSSION 
observed in the tail or null region. For a cold gas coming in at 100° or something like 
that, you find that the absorption cross section at the center of the line is approximately 
100 times larger than the average absorption cross section calculated in the optically thin 
approximation. In general this makes a very small difference in all calculations with 
regard to Lyman a. and so the optically thin approximation is very good. But this differ- 
ence has two effects. The first is that it can produce a variable addition to the gravita- 
tional field from the sun so that the field may change from inward at larger distances to 
an outward field at smaller distances. This would effectively give a large temperature 
effect and throw some particles into the null region and thus give you a small Lyman a 
flux there. I don’t know the magnitude of this effect, so I don’t know whether it’s 
enough to account for the observations. But from the calculations I made with regard to 
something else I think it may be. 

My second point is also with regard to the optically thin approximation for the 
absorption at the center of the line. Even though the gas is generally optically thin, this 
absorption will tend to increase the density that you would predict. At 0.1 particle/ cm 3 
the increase can be as much as a factor of 2. Consequently, instead of from 0.03 to 
0.12 cm" 3 , your density might be from 0.06 to 0.24 cm" 3 , As you may have noted, this 
difference makes a considerable difference in the' effect of the interstellar gas on the solar 
wind. 

G. E. Thomas In regard to the first point concerning the effect of absorption on the 
solar Lyman a as it progresses outward, I would point out that this absorption is due to a 
scattering by hydrogen atoms so that eventually the center of the line is eaten out by 
neutral hydrogen. But, as I recall, the mean free path for Lyman a proton with these sorts 
of densities is of the order of 10 AU or so. I can’t remember exactly . We decided that 
such an effect was really unimportant, that is, out to about 5 AU or so. Such things might 
happen beyond 5 AU, but they would not strongly affect the Lyman a scattering mea- 
surements. 

Now let me comment on the second point concerning the importance of multiple 
scattering. If you again take the mean free path to be about 10 AU, it turns out that 
multiple scattering effects produce an intensity that is small compared to the singly 
scattered intensity we calculated. If you turn the argument around and suppose there’s a 
lot of multiple scattering! then you quickly come to the conclusion that you shouldn’t 
get the large asymmetries that are observed. As you see the ratio between the maximum 
and the minimum is better than 2/1. If multiple scattering were dominant, it would 
tend to wipe out that difference between the maximum and the minimum. 

T. E. Holier I was talking about absorption near the center of the line with regard to 
the second point as well as the first. With regard to the second point, it does make a 
factor of 2 difference. With regard to the first point, it may or may not be significant; it’s 
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hard to say. 

G. E. Thomas Those are difficult calculations. 

J. C Brandt What is the absorption cross section in the center of the line? 

T. 2s. Holzer I can’t give you a number, but if you assume that the Lyman a line has 
a flat topped peak 1 A wide and then you assume the temperature of the neutral gas to be 
100° you can calculate the cross section of the center of the line and you can calculate a 
mean cross section for the whole line. The cross section of the center of the line is a 
factor of 100 larger than the mean cross section. Actually, it’s a factor of 300 larger for a 
temperature of 100° and the other assumptions. 

J. C. Brandt Take it at the center line, and you’ll get the maximum value. 

T. 2s. Holzer Right. But the point is that there is a factor of 100 or 300 between the 
average cross section and the central cross section, because in taking the mean cross 
section one tends to neglect the effects of any absorption of the line. In any case, I 
wanted to suggest that the actual absorption effects are somewhat important because, 
although they don’t change Dr. Thomas’ results basically, they can increase the density 
by a factor of 2. 

G. E \ Thomas I think the importance of multiple scattering depends very much on 
the temperature. For the hot model, the cross section at the center of the line is much 
smaller than in the cold model, and mean free paths are larger so the net effect is to 
reduce the importance of multiple scattering. For reasons I have already given, I tend to 
believe the hot model. 

F. Scherb I would like to mention a likely perturbation on the picture that you 
presented. As I understand it, the axis along which the flow is directed lies nearly in the 
ecliptic plane. If I have done the arithmetic correctly, about a year from now Jupiter will 
cross this axis at 5 AU; thus, Jupiter, a source at 3 AU, and the sun will lie on a line, and 
the gravitational perturbations by Jupiter might be observable if one were observing the 
Lyman a distribution in this region of the sky. The perturbations might not be entirely 
negligible in view of the fact that the sun’s gravity appears to be virtually absent as far as 
the hydrogen is concerned. 

G. E. Thomas Well, it might be interesting to look for a bright spot and I agree, the 
bright spot might be there, but I’m not sure what we would learn by it. 

J. C. Brandt Jupiter’s shadow also might be quite long. I’m sure it couldn’t be done 
with the wide-resolution instruments on OGO 5. Instruments with better spatial resolu- 
tion might see the shadow. I’m sure it’s quite long, approximately 1 AU. 

C. P. Sonett Just one very short note, which I think bears on this. There is a 
considerable body of work now indicating that the temperature of some interstellar 
clouds is quite high, perhaps due to heating by low-energy cosmic rays. 

J. C. Brandt You can get a good argument on that particular point with any group 
of astronomers containing two or more people. Some of us have been discussing a process 
of heating the interstellar medium with 304A (He//) photons from a supernova event. 
The calculations indicate that the process contributes. However, in all of the colloquia 
and discussions I’ve heard on this subject recently, people just throw up their hands at the 
end. Ultimately, this will be sorted out, but I don’t think there’s any agreement at the 
present time. 

W. I. Axford Could I ask if anyone knows what is the status of observations of the 
temperature of the neutral hydrogen, apart from the 21 -cm measurements? Is there 
anything based on the width of absorption lines that suggests there might be a high 
temperature somewhere? 

/. C. Brandt If you look, for example, at the spectra, of the kind that Munch takes, 
or used to take, the lines are very narrow and very closely spaced and in fact very hard to 
tell apart. The gas appears to be cool. 



E. N. Parker The only comment I could add is that if one sees a broad absorption 
line, he attributes it to turbulence rather than to the high temperature of the neutral 
hydrogen. So I don’t know that one can unravel it by standard observations apart from 
those of the 21 -cm lines. 

/. C. Brandt Yes, there is the random component of 1 0 or 20 km/sec. 

W. C. Feldman Concerning what you would learn by looking at Jupiter when it 
comes around, I would like to point out that you have the same sort of gravitational 
focusing effect as was shown by Axford, and the only thing that limits the cusp is 
temperature. By looking at the size and brightness of the spot you might get an indication 
of the temperature of the gas. 

F Scherb With regard to the temperature of the neutral gas there is a program under 
way now at Wisconsin to measure Balmer emission lines from the neutral gas due to 
recombination of ionized hydrogen. By making measurements of H a and you can in 
fact infer the temperature of the gas because the probability of emission of an photon 

is virtually independent of the temperature at which the recombination occurs, whereas 
this is not true for H a . This effect is most useful for temperatures below 1000° K. Above 
1000° one just measures the line profile. So it should be possible, with the proper 
sensitivity and spectral resolution, to infer the temperatures by studying the H a and 
emission from the clouds. 

J. C Brandt Where in the galaxy do you expect to get sufficient brightness to be 
able to observe this? 

F Scherb This has already been observed by looking just off the Crab, far enough so 
that the Crab itself was not in the field of view. The Hj3 emission was detected. The 
program will be continued, with better equipment, at the Goddard Space Flight Center 
over the next 2 years. 

W. L Axford One way of getting an answer to this question of the temperature 
would be to look at the scattered radiation from helium atoms at 584A. Helium atoms 
are not affected by radiation pressure and their thermal speed is half that of hydrogen 
atoms. Therefore, if the gas is relatively cold, the focusing effect will be quite pro- 
nounced. Although one expects that the general intensity from the helium will only be 
perhaps 1 to 10/? depending on the width of the solar line, the bright spot associated 
with the helium cusp could be very much more intense if the gas is cold. It might well be 
100/? or so in a small area, and I believe that this could be detected quite easily. 



